Nonlinear phenomenology from quantum mechanics: Soliton in a lattice 



Juha Javanainen and Uttam Shrestha 



o 
o 



(N 



43 

9 L,: 



> 

On 

o 

cn 
o 
o> 
o 



X 



Department of Physics, University of Connecticut, Storrs, 
(Dated: March 29, 2009) 



CT 06269-3046 



We study a soliton in an optical lattice holding bosonic atoms quantum mechanically using both 
an exact numerical solution and Quantum Monte Carlo (QMC) simulations. The computation of 
the state is combined with an explicit account of the measurements of the numbers of the atoms at 
the lattice sites. In particular, importance sampling in the QMC method arguably produces faithful 
simulations of the outcomes of individual experiments. Even though the quantum state is invariant 
under lattice translations, an experiment may show a noisy version of the localized classical soliton. 
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Quantum mechanics is widely accepted as the funda- 
mental framework for physics. Nonetheless, a variety of 
nonlinear theories have met with unquestionable success 
in the modeling of macroscopic (many-body) systems. 
This state of the affairs begs for the question of how linear 
quantum mechanics mimics nonlinear behavior. Consider 
as an example a soliton of light propagating in a nonlinear 
optical fiber. Full quantum theory has solutions that in 
many respects resemble classical solitons but finding 
such a quantum soliton does not solve the problem of lin- 
ear versus nonlinear dynamic. For instance, in a transla- 
tionally invariant optical fiber the quantum ground state 
may be chosen to be translationally invariant. Why is 
it, then, that localized solitons are so strongly favored by 
Nature that they form seemingly spontaneously? 

It has been suggested that measurements play a part in 
the answer 0], but our tenet goes much further: Opera- 
tionally, nonlinear phenomenology is created by the pro- 
cess of measurement. In this Letter we demonstrate our 
views by discussing bosonic atoms in an optical lattice, 
a problem that commands much current interest in its 
own rig ht 0,0] and also presents a discrete-space version 
of the soliton in nonlinear optics. We solve the quantum 
mechanical ground state either by exact diagonalization 
or by Quantum Monte Carlo (QMC) simulations [US], 
but with a measurement theoretical twist. Namely, while 
a run of a QMC simulation produces the numbers of the 
atoms at each lattice site to be used in calculations of 
various ground-state expectation values, we argue that 
these occupation numbers are also a faithful simulation of 
what one would find in a single experiment that measures 
the occupation numbers. When classically one expects a 
soliton, each run of the QMC simulation, and hence each 
individual experiment, also shows a distribution of the 
atoms over the lattice sites like a classical soliton. 

Quantum mechanically, we have the Bose-Hubbard 
model with the Hamiltonian 



~^L+A + bl-ibk) + 'i'blblhbk 



(1) 



is the amplitude for site-to-nearest-site tunneling, and k 
characterizes on-site atom-atom interactions. Our L-site 
lattice is periodic so that k — L and k — are the same 
site. This boundary condition, which would physically 
correspond to a ring lattice, is specifically chosen to make 
the Hamiltonian invariant under lattice translations. The 
total number of atoms N — J^k b k b k is a constant of the 
motion, and its value is denoted by N. 

The corresponding classical Hamiltonian is 



■^(b* k+1 b k + bl^bk) 



fN 4 



(2) 



Here bk are complex numbers such that \b k \ 2 stands for 
the number of atoms at the site k. The quantities bk and 
bt> are regarded as classical canonical conjugates with 



the Poisson brackets {bk,b 
equation of motion for bk 

S 



-(i/h)Sk,k'- Hamilton's 



ibk = ~{bk+x + b k -i) + K\b k \ 2 b k , 



(3) 



Here bk annihilates an atom at the site k = 0, . 
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the lattice analog of the Gross-Pitaevskii equation, is 
commonly called Discrete Nonlinear Schrodinger equa- 
tion (DNLSE). The total atom number N = J2 k M 2 is 
again a constant of the motion. 

The ground state of the classical Hamiltonian |2} is 
found, e.g., by integrating the DNLSE © in imagi- 
nary time. For repulsive interactions between the atoms, 
k > 0, the lowest-energy state is spatially uniform, 
bk = VNe lip , where ip is a global phase. On the 
other hand, for a sufficiently strong attractive interac- 
tion k < 0, the lowest-energy state is a soliton localized 
around some lattice site Q . A soliton is a stationary state 
of the DNLSE; the amplitudes evolve in time according 
to b k (t) = e- Hlt / h b k (0), where ji is the chemical poten- 
tial. Physically, the soliton reflects a balance between 
the tendencies of the atoms to collect together drawn by 
the attractive interactions, and to disperse as a result 
of the site-to-site hopping. A lattice-translated solitonic 
ground state is also a ground state, but a nontrivial su- 
perposition of ground-state solitons is not. Nonlinearity 



2 



makes localization an inescapable feature of the solitonic 
ground state. 

In search of the quantum counterpart of a soliton 
we study the Hamiltonian ([1]). Introduce the ther- 
mal density operator at the inverse temperature /?, 
p(j3) = e~ l3H /Tr[e~P H ], the number operators of the 
sites k, hk = bVbk, the collection of these operators 
n = {^o» • • • ) fi/L-i}, and the simultaneous eigenstates 
|n) of the number operators fi characterized by a collec- 
tion of L integers n = {no, . . . , ri£_i}. We call such a set 
n of the occupations numbers n^ a number state of the 
lattice. As the operators hk commute, according to stan- 
dard quantum mechanics their values can be measured 
simultaneously. We say that the operator n is measured. 
Any function /(n) of the number state defines in a nat- 
ural way the operator function /(n), and the thermal 
expectation value of /(n) equals 

(/(n))=5>(n)/(n); P(n) = J^>L . (4 ) 

On the other hand, suppose that by some means, call it 
importance sampling, it is possible to produce number 
states n m , m = 1,2,... at random, but with frequencies 
proportional to the corresponding P(n m ) as in Eq. ([J|. 
The expectation value in Eq. (T5)) may then also be written 

M 

</(n)> = </(n)> = lim - £ /(n m ) . (5) 

m— 1 

Considering any quantity that can be expressed in 
terms of expectation values of functions of the form /(n), 
the values from quantum measurements and from impor- 
tance sampling agree. This applies to average atom num- 
ber at a site, correlations of atom numbers between the 
sites, and so on: The measured values of the number op- 
erators hk agree in every possible statistical characteriza- 
tion with the occupation numbers in the importance- 
sampled number states n m . One experiment and one in- 
stance of sampling each produce a random set of occupa- 
tion numbers, but in repeated runs the statistics of these 
random numbers could be verified to be the same for the 
experiments and for importance sampling. We therefore 
formulate our key interpretative hypothesis: Each number 
state n m coming from importance sampling is a represen- 
tative outcome of an experiment measuring the operator 
h. 

For brevity we focus on the low-temperature limit, 
(3 — > oo. If the ground state of the Hamiltonian (TTJ) is 
nondegenerate, the density operator is simply the pro- 
jection onto the ground state. On the other hand, if the 
ground state is degenerate, the zero-temperature state is 
an equal mixture of the degenerate ground states. As a 
lattice translation is a symmetry operation of the Hamil- 
tonian, it is possible to select all energy eigenstates so 
that they are invariant under lattice translations. Hence, 



0.3 



Si 0.2 



0.1 



0.0 







+ 


N =16 






+ 


+ 




64 


4- 


+ 




A 


+ 


256 + 


A 




+ 


/ 


\ 


1024 / 

/ 




+ 






\ + 


+ / . 




\ 

\ 


/ 




! \ 


/ ; 




/ 

/ 

/ 




i \ 

! \ + 
\ 


/ ; 
+ / ; 
+ / 


\ \ 


'; \ 
■; \ 


+ / f 

1 ; 


1 


'', \ 

\ 

\ 


/ 

/ ! 

/ 










\ 














/ 










\ V 








JJl 

























0.0 0.2 0.4 0.6 0.8 1.0 



n/N 

FIG. 1: Representation of the probabilities P n that n atoms 
are in one given site of a two-site lattice for different numbers 
of the atoms N , as indicated in the legend. This figure is from 
exact quantum solutions for the fixed value of the parameter 
Nk/5 = -2.309. 

the unique zero-temperature density operator must be 
invariant under lattice translations. Unlike the classical 
soliton, the quantum mechanical ground state (or, for 
that matter, any thermal-equilibrium state) will not sin- 
gle out any particular lattice site. 

We first take up the two-site lattice, L = 2. Simple as 
the problem is, for parameters such as Nk/5 = —2.309 
we already have the conundrum that the classical ground 
state is a soliton with | of the atoms in one site and \ 
in the other, whereas the quantum ground state is trans- 
lation invariant. So, what would the experiments see? 

Here the quantum mechanical state space is spanned 
by the N + 1 vectors \n,N — n), where n (N — n) is the 
number of the atoms in the site k = (k = 1). The two- 
site problem is trivial to solve numerically. For a fixed 
value of the parameter Nk/5 that classically gives a soli- 
ton, and in the limit N — > oo, the ground state is dou- 
bly degenerate. Therefore we write the low-temperature 
ground state as a 50/50 statistical mixture of the two 
lowest-energy quantum states. As far as observed atom 
numbers are concerned, the complete statistics is deter- 
mined by the probabilities P n that n of the atoms are 
found in the k = site. 

In Fig. [1] we plot the quantity \f~N P n as a function of 
the fraction of the atoms in the k = site, n/N . Here 
Nk/5 = —2.309 is fixed, and the atom number N is 
varied as indicated in the legend. Importance sampling 
would produce number states n = {n, N — n} at fre- 
quencies proportional to P n . As Fig. [1] shows, with an 
increasing atom number N the sampling, and the exper- 
iments, would give an increasingly accurate split of the 
atoms so that | of the them is in one site and j in the 
other, just like in the classical soliton. The reason for the 
y~N multiplier for the probabilities is that one may then 
easily see from Fig. [T] that the fluctuations in the 3 : 1 
split scales with the atom number N like 1/yN. 

While the Hamiltonian is invariant under a lattice 
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translation, an individual measurement result of the 
atoms numbers with the approximate 3 : 1 split between 
the sites is not. However, either one of the two sites 
is the one with the higher occupation number with the 
same probability, so that averaging over the experiments 
restores the symmetry. 

The resemblance of this situation to spontaneously 
broken translation symmetry has been brought up in nu- 
merous papers on the continuous soliton problem and 
the attendant connection to measurement theory has also 
been noted but here we take a very direct approach. 

Consider interference between two Bose-Einstein con- 
densates, an analogous problem from the past [1, 0, [1(3] 
that continues to evoke new angles The traditional 
view is that gauge symmetry of the condensate is broken, 
which endows the BEC with a phase. The difference in 
the phases is observable when two condensates are made 
to overlap, which indeed produces an interference pat- 
tern [12 ]. In contrast, we have predicted an interference 
pattern for two number-state condensates without ever 
assuming any broken symmetry or phase [§]. The key 
was to simulate the measurements of the positions of 
the atoms. There are correlations between atomic po- 
sitions in the state vector of the two overlapping conden- 
sates. Every time the position of one atom is observed, 
the state gets reduced so that it is compatible with the 
newly gained measurement result. This reduction brings 
out the correlations by modifying the position distribu- 
tion for the next atom to be observed. Continuing in this 
way, the measured positions of the atoms produce an in- 
terference pattern. It comes about as a combination of 
two elements: correlations between atomic positions in 
the state, and measurements that convert the correla- 
tions into the observations. 

Easy as the corresponding measurement simulation 
would be in the case of a two-site lattice 13, 13] , we have 
not carried it out because the analog is clear enough as 
it is. Even if the ground state is translationally invari- 
ant, the potential for a soliton is there in the correlations 
between the positions of the atoms: The attractive inter- 
actions favor atoms collecting in the same site. However, 
it is the observations that ultimately make the soliton. 

With increasing numbers of atoms and lattice sites an 
exact numerical solution 14j of the lattice system be- 
comes impractical. Besides, the solution does not tell 
what an experiment would see; one has to add a measure- 
ment simulation along the lines of Refs. [3, llfj, ll3j • QMC 
simulations may provide a sweeping answer to the ques- 
tions of both finding the state and measuring it. This 
is because importance sampling is a core idea in QMC 
methods, and in Monte-Carlo simulations in general 15 1. 

Our QMC simulations take place on a L x Np grid 16| , 
where the first dimension represents the lattice sites and 
the second dimension corresponds to Np steps in inverse 
temperature with the size f3/Np. We use an elementary 
world line algorithm combined with the checkerboard de- 




FIG. 2: Atom numbers at lattice sites k for the classi- 
cal soliton (solid circles) and for five solitons sampled from 
quantum mechanics; the symbols for the six solitons are all 
resolved at the site k — 9. The parameters are L = 16, 
N = 256, and k/S = -0.004. The temperature is 5/100 in 
units of frequency. Given the occupation numbers for the soli- 
ton 71% and for the QMC results , the £ 2 distances between 
the soliton and the five simulations, d = Et;( n fc — n k^ 2 ^^ 2 
equal 12.8, 20.9, 18.8, 17.8, and 15.8. 



composition of the Hamiltonian [3, Q. The grid is ini- 
tially seeded with the same number state of the lattice at 
each step of the inverse temperature. A site on the grid is 
selected in a suitable way, say, at random, and at the site 
a move is proposed, also at random; an atom moves left 
or right to or from an adjacent lattice site. The move is 
either accepted or rejected with a probability that is cal- 
culated from the state of the grid in the immediate neigh- 
borhood of the chosen site. Moves are attempted in this 
way until the grid has relaxed to a steady state and ex- 
ecutes fluctuations around it. We use periodic boundary 
conditions in the inverse-temperature direction as well, 
so that a number state n in the relaxed grid at any step 
of the inverse temperature qualifies as a sample drawn 
at a frequency proportional to P(n), Eq. ((4|). Overall, a 
QMC simulation sets up an artificial dynamics that has 
nothing to do with the real dynamics of the lattice, but 
realizes importance sampling for the thermal state. 

To demonstrate spontaneous breaking of translation 
symmetry we would like to initialize the grid so that 
every site has the same atom number. This works for 
small numbers of atoms and lattice sites, and for suitable 
strengths of the attractive interaction solitons are found 
once the simulation has converged. However, a soliton 
is nucleated from fluctuations of atom numbers in the 
grid, and for large atom numbers a sufficiently large fluc- 
tuation does not occur during the number of attempted 
moves we are willing to wait out. In such a case we form 
the soliton by seeding the grid in the direction of the 
lattice sites with a narrow distribution of the atoms. 

Let us now compare the classical soliton and results 
from quantum simulations for L = 16 lattice sites and 
N = 256 atoms. We choose the interaction strength 
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k/5 = —0.004, for which there is a classical soliton that 
spans several sites but is well contained inside the lat- 
tice length L. In Fig. Owe plot both the classical soliton 
(solid circles) and five importance-sampled solitons from 
the QMC method (various symbols, best resolved at the 
lattice site k = 9). QMC results have been translated 
along the lattice so that the root-mean square difference 
from the classical soliton is minimized. The fluctuations 
in the simulations are rather large, but an eye would rec- 
ognize every result of sampling as a soliton. 

We have discussed explicitly two examples summarized 
in Figs. Q] and [2j but in our experience the results are 
generic: If for some fixed values of L and Nk/5 there is 
a classical soliton, in the limit of large N a measurement 
on the quantum system will produce a soliton as well. 

Our main technical point that each outcome from im- 
portance sampling faithfully represents an outcome of a 
single experiment may seem self-evident, but in practice 
it appears to be difficult to grasp this type of an argu- 
ment. We therefore expand on the underlying philosophy. 

Ultimately, the claim is that the statistics of a random 
variable, here n, determines how the realizations look 
like. 



17, 1 



We have resorted to the same idea before 
for instance when we predicted "quantum jumps" in the 
intensity of the light scattered off a three-level system 
with a long-lived "shelving state" ab initio from quan- 



tum mechanics 17] . The intensity correlation functions 
for the scattered light are continuous functions of their 
time arguments, but they are the correlation functions 
of a Markov jump process. Hence, the light intensity 
should look like a Markov jump process, i.e., execute 
abrupt jumps. These quantum jumps, originally envis- 
aged by Dehmelt [l9j], were subsequently observed [2(| 
and are now a standard laboratory technique for detect- 
ing weak transitions. In the present case we have no prior 
experience with the underlying random variable, but im- 
portance sampling shows that the realizations look like 
solitons. Hence, experiments should produce solitons. 

To further illustrate the paradigm shift in our reason- 
ing we compare with traditional QMC simulations. Or- 
dinarily one would produce a number of samples and use 
them to calculate average quantities such as occupation 
numbers of the sites, standard deviations of the occu- 
pation numbers, correlations in atom numbers between 
adjacent sites, and so on; and maybe compare such av- 
erages with averages extracted from experiments. How- 
ever, such a process may miss the forest from the trees: 
Even though the results from importance sampling dis- 
play solitons, solitons may be difficult to uncover from 
the averaged quantities. 

Recognizing the soliton as the basic outcome from mea- 
surements should inspire novel theoretical work. Con- 
sider analyzing fluctuations of a soliton using lineariza- 
tion of quantum mechanics around the classical soliton 
2l| . which in the context of Bose-Einstein condensates 
is commonly called Bogoliubov theory. For a fixed num- 



ber of sites the DNLSE can be rewritten so that the 
only dimensionless parameter containing the atom num- 
ber is Nk/6, the same applies to the linearization of the 
DNLSE that gives the Bogoliubov theory, and therefore 
also to atom number fluctuations around the classical 
soliton [2^, . We have here the interesting situation 
that the atom number fluctuations in the Bogoliubov the- 
ory can only depend on atom number through the param- 
eter Nk/6, which directly contradicts our Fig. [TJ 

There are many layers in our discussion. At the lowest 
level we introduce the measurement-theoretical idea that 
importance sampling in our simple QMC method liter- 
ally simulates the experiments. New ways of comparing 
theory and experiment then open up. We apply this line 
of thought to predict localized solitons for atoms in an 
optical lattice from quantum mechanics in spite of the 
fact that the quantum state is translationally invariant 
and does not favor any particular lattice site. At the 
top level we promote the notion that classical nonlinear 
phenomena in a macroscopic system are implicit in the 
correlations within the quantum state, but measurements 
are the agent that ultimately brings them forth. 

This work is supported in part by NSF (PHY- 
0750668). 
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